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Abstract 

The Dirac equation is considered in the background of potentials of several 
types, namely scalar and vector-potentials as well as "Dirac-oscillator" potential 
or some of its gener alisations. We investigate the radial Dirac equation within 
a quite general spherically symmetric form for these potentials and we analyse 
some exactly and quasi exactly solvable properties of the underlying matricial 
linear operators. 







1 Introduction 



In quantum mechanics, the Schròdinger equations which can be completely solved 
by algebraic methods are rather exceptional. One can attenuate the condition of 
complete solvability by asking that, at least, a few eigenvectors can be obtained 
by an algebraic method. The so called "quasi exactly solvable" (QES) equations 
[1, 2] refer to a class of quantum hamiltonians which possess precisely this property. 
The corresponding operators can be set in correspondance with finite dimensional 
representations of some Lie algebras and, accordingly, a few of their eigenvectors can 
be computed by solving an algebraic equation. 

It turns out that QES Schròdinger equations often occur as suitable extensions of 
exactly solvable ones. The most famous example [1] is the one dimensional quantum 
harmònic oscillator. When completed by a suitable quartic plus sextic potential the 
equation of the harmònic oscillator becomes QES. In three dimensions, the proto- 
type of exactly solvable quantum hamiltonian corresponds to the charged particle 
in a central Coulomb potential. When the Coulomb potential is supplemented by a 
confming potential of the form V = u 2 r 2 + Àr, the equation becomes QES [3]; the 
corresponding spectral problem is, however, of type II [1]; that is to say that the 
energy eigenvalue of the initial Schròdinger problem is not the spectral parameter 
of the algebraic equation. A specific combination of the physical coupling constants 
plays the role of the spectral parameter of the QES equation; one energy level can be 
determined for the corresponding vàlues of the coupling constants. 

The possibility of obtaining operators with interesting algebraic properties can 
also be investigated for relativistic equations like the Dirac equation which leads, 
in general, to a system of coupled equations. Although the classification of QES 
equations is well understood for scalar equations, this is not the case for systems and 
a complete classification of matrix-valued QES operators is still missing. As a general 
rule the criteria to fullfill the QES property are more severe for systems of equations 
and it turns out interesting to consider "type II" systems as well [4]. The family 
of QES problems obtained in [4] was extended in [5] whom proposed a method to 
incorporate screened Coulomb and scalar potentials. 

More recently the Dirac-Pauli equation was reconsidered [6] in the context of 
neutral particle interacting with an electromagnètic field and shown to be strongly 
related to the case of a Dirac oscillator. A similar problem was emphasized is 1+2 
dimensions in [7] . With suitable form of the radial potentials, the corresponding Dirac 
equation can be quasi exactly solvable. 

In this paper, we investigate the occurrence of explicit solutions in the framework 
of the Dirac equation coupled to a class of external radial potentials, extending the 
cases of Dirac oscillator and the Dirac-Coulomb problem and generalizing the choices 
of [6, 7]. The general physical problem to deal with is presented in Sect.2. Three 
cases for which the underlying operator is completely solvable, namely the Dirac 
oscillator, the extended Dirac oscillator and the Dirac-Coulomb problem are presented 
in Sect.3; the emphasis is put on the construction of equivalent operators preserving 
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an infinite flag of vector spaces of polynomials. In both cases the quantization of the 
eigenenergy appears as a necessary condition for these invariant spaces to exist. The 
general problem of constructing QES operator out of the generalized Dirac oscillator 
is adressed in Sect. 4 

2 Dirac equation and radial potentials 

We will study the radial equations associated with the 3+1 dimensional Dirac equa- 
tion coupled to (radial) scalar and vector potentials. We start with the Dirac-Pauli 
equation : 

{{id, - eAJ>f - i/in^i^ — M — W(r))^(x ,x k ) = (1) 

where is a vector potential, W is a scalar potential and F^ is the electromagnètic 
field. The charge, mass and anomalous magnètic moment of the spin-1/2 particle are 
respectively denoted e, M, jj, n . 

After standard manipulations, namely the separation of the time variable and 
the separation of the angular variables in a central potential, the radial equation 
associated with (1) takes a conventional form of a 2*2 matrix equation [6]. Here we 
will assume the quite general form 

i-* + Hn(B-E) M-e-V + W )(f) =Q m 
M + e + V + W ± + * + l· i n (B + E) ) \g ) [) 

where e is the energy parameter while k is the total angular momentum. In this 
equation we have set four independent external radial potentials. The parts W(r) 
and V(r) denote respectively a scalar potential and the time part of vector potential 
Ap. The part E(r) and B(r) are related respectively to a radial elèctric and magnètic 
fields which couple through the anomalous momentum fj, n of the particle. In the 
following we will consider the mathematical problem of finding explicit solutions of 
the above equation with the form of the generalized potentials E(r), B(r), V(r), W(r) 
fixed a priori by hand. This general point of vue makes sense since e.g. the Dirac 

— * 

oscillator is recovered for fj, n E — > mur. 

In the following, we will restrict ourselves to radial potentials of the forms 

a s 6 s s 

V(r) = - + 5>ir» , W(r) = ^ + Y / frr t , E(r) =^7^' (3) 

r i=l r i=l i=0 

along with [6] we assume B(r) =0. 

3 Exactly solvable cases 

In this section we construct explicitely the infinite flag of invariant vector spaces for 
two particular cases for which the operator above turns out to be exactly solvable. 
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The technique is such that the quantization of the energy comes out as a necessary 
condition for the operator to preserve a vector space of the form V(n) © V(n — 1), 
where V(k),k G N denotes the polynomial of degree ot most k in an appropriate 
variable, say x. 

Because the vector space V(n) © V(n — 1) constitutes the basis of a particular 
representation of the graded- àlgebra osp(2,2) [8], these results demonstrate that the 
operator H corresponding to these cases is indeed equivalent to an element of the 
envelopping àlgebra of osp(2,2). 

3.1 Dirac oscillator 

For this case one has V = 0, W = 0, E(r) = —r. The invariant spaces of the 
corresponding operator H are revealed in terms of spaces of polynomials after we 
perform the following transformation 

H = A- l HAC (4) 

with the operator A and constant matrix C defined according to 

Choosing 9 = k and u 2 = ji n and using x = r 2 as a new variable, the Hamiltonian 
takes the form 

H = x( ° °)+2x±( 1 M ~ £ ) (6) 

X { M + e ) +ZX dx { 1 ) + \0 2k-i) W 

which turns out to be a linear combination of the generators of the súper àlgebra 
osp(2,2) (in the suitable representation first pointed out in [8]), provided the energy 
parameter e is of the form 

e 2 = M 2 + 4n/i n , n G N (7) 

reproducing the spectrum of the Dirac oscillator equation (see e.g. [12]). Indeed, in 
this case, the matrix element H í2 becomes proportional to the operator x-^ — n and 
correspondingly, the operator preserves the vector space of polynomials of the form 

(p n _i(a;), 

This way of obtaining the spectrum of the Dirac operator by enforcing it to pre- 
serve an infinite family of (finite dimensional) vector spaces reveals the hidden alge- 
braic structure of the Dirac oscillator equation by its relation with osp(2,2). 

3.2 Extended Dirac oscillator 

The result of the Standard Dirac oscillator discussed above can be extended to a case 
of type II exactly solvable system if the various potentials are suitably chosen. In the 
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purpose to exhibit a such possible extension, let us consider the potentials with the 
following form : 



V(r) = - , W(r) = (3 ir + ^ , E(r) = 7o + 7i r . (8) 
r r 

In order to reveal the hidden àlgebra structure of the corresponding H several changes 
of variables and/or functions are necessary. First, it is convenient to write 71 = 
_Rcos(2u;), (3i = i?sin(2a;).Secondly, we "gauge rotate" the operator by means of 

U Sexp(-Rr 2 /2-Tr) , S = diag(rV) , U = ( C ° SUJ ~ sinu; ) (9) 
v ' ' y ' y sinw cosa; J v ' 

Then, the following vàlues have to be imposed (with c = cos(2u;) , s = sin(2a;)) 

a = , (3 = -Ktan2u , 7o = -Mtan2cu ,T = , 6 = kc , 6-6 = A = -1 (10) 

and if the so obtained operator, say H, is further transformed by means of 

and the variable x = r 2 is used, we finally obtain an operator which preserves 
(j) n -i(x) , Pn( x )Y ■ The energy quantization formula reads in this case 

M 2 

E 2 = —— + 2nR , )i6R (12) 
cr 

Finally, the quantum Hamiltonian obtained in this way generalises the conventional 
Dirac oscillator by means of the addition of the angle uj appearing through the cou- 
pling constants /3i,7i- All the other coupling constants (namely 7o,/3,a) are fixed in 
terms of M, k and u. Since a = turns out to be a necessary condition, the attempt 
above further reveals the difficulty (if not the impossibility...) to obtain a hidden alge- 
braic form of the combined relativistic problem Dirac Oscillator + Coulomb potential. 
In the following section, we will succeed in producing partial algebraic solutions to 
this problem. 



3.3 Dirac-Coulomb problem 

This case is well known and was presented in [4] adopting the point of view of exactly 
solvable operators. However we present briefly the construction here for the purpose 
of completeness and for the sake of comparaison with the case studied in the previous 
section. The conditions on the potentials are — — 0, E — B — 0. The relevant 
transformation of the Hamiltonian reads 

H = AU- l HUC (13) 
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U = /exp(-Ar) f U ~ U ~ , u ± = VM±e (14) 

V u + ~ u + I 





with 



(15) 

Choosing the arbitrary parameters appearing in the operator U according to 

2 = « 2 + /3 2 -a 2 , X 2 = M 2 -e 2 (16) 

and multiplying the equations by r, leads to the new operator H which reads (up to 
diagonal elements which depend on r J;) 





(17) 



with 



^vm^-^^ 2 - 02 (18) 



If the quantity i] is imposed to be an integer, n, the operator H manifestly pre- 
serves the vector space P(n — 1) © P(n) for n — 1, 2, 3, . . . (as before P(n) denotes 
the space of polynomials of degree less or equal to nin r). 

The condition rj = n leads to the quantization of the energy of the system [9, 10]. 
The celebrated spectrum of the relativistic hydrogen atom (in the case /3 — 0). As for 
the case of the Dirac oscillator it is obtained by requiring that the reduced spherically 
symmetric Hamiltonian can be written in terms of the generators [8] of the súper 
àlgebra osp(2,2). The infinite series of finite dimensional vector spaces preserved by 
this realization corresponds to the eigenvector of the problem. 



4 Quasi exactly solvable cases 
4.1 Planar case 

Here we consider the case of a planar Dirac electron in Coulomb plus magnètic field. 
This case was studied in [7, 11]. In the general framework of Eq.(2) it corresponds 
to B(r) = 0, W(r) = and E(r) — > eBr/2. We also note two main differences that 
r should now be interpreted as a two-dimensional radial variable and that k now 
takes half-integer vàlues. Also, in the following discussion we will rescale the radial 
variable according to x — rVeB. Posing, along with [7], / = x' 1 exp(—x 2 /4)Q(x) , 
g = x~ t exp(—x 2 /4:)P(x) where P, Q are polynomials leading to the following equations 
for P, Q : 

P' + ^P-(e-M+^)Q =0 , (19) 
Q' + ^Q-xQ + (e + M + ^)P =0 (20) 
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Where e and M have been rescaled appropriately. The polynomials Q, P should have 
degrees n, n + 1 respectively, nevertheless the number of equations 2n + 5 exceeds by 
one unit the number of free parameters 2n + 4 and therefore a algebraic solution only 
exist for specific vàlues of a. Two of the equations allow to fix immediately the vàlues 
of 7 and e according to 

7 = V/t 2 - a 2 , e 2 = M 2 + 7 + /t + n + l (21) 

Unfortunately, we did not manage to find the extra condition for a in a simple way 
for genèric vàlues of k, M and n. 

In the simplest case, n = 0, we find 

e = - 1 -(M±VWT2) , a 2 = -I±|^ (22) 

showing that two branches of vàlues of a lead to polynomial solutions. However 
these eigenvectors are available only for k < 0. In particular, this indicates that the 
cylindrical symmetric ground state is not one of the algebraic solutions. 

Already forn = 1 the final equation relating a, M, k is very involved. The allowed 
vàlues of a with M, n fixed can be determined numerically. For the case k — 1/2, 
there are four possible vàlues of a and two corresponding vàlues of the energy e. They 
are presented on Fig.l for Mg [0,1]. For k — 3/2 the picture is qualitatively similar 
to Fig.l. 

4.2 Extended Dirac oscillators 

The case treated in the previous section is physically important, but we have seen 
that no algebraic solution can be constructed for genèric vàlues of the Coulomb and 
oscillator coupling constants, respectively a, B. We would like to construct a model 
where algebraic solutions exist for genèric vàlues of the Coulomb and oscillator con- 
stants. In this purpose, we consider Eq.(2) with an extended choice for oscillator's 
parameters. Namely, we set 

OL 

V(r) = - , W(r) = far E{r) = 7o + 7i r . (23) 
r 

This choice constitutes one possible relativistic generalisation of the non-relativistic 
problem considered in [3] Here the Coulomb interaction is supplemented by two types 
of confining (i.e. linear in r) interactions. Note that, adding a constant to the 
potential W(r) is equivalent to redefine the mass M. To the contrary, the constant 
7o cannot be eliminated by a redefinition of the physical quantities. 
We look for solutions of the form 

(^)=r e exp(-^r 2 -A 1 r)^) (24) 
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where 9, Ai, À2 are constants and p(r), q(r) are chosen as polynomials in r. Inserting 
(24) into (2), we obtain the counterpart of (17) with 

_ í D - X 2 r 2 - A x r + 9 - k + 7 r + ^r 2 ^r 2 -a + r(M - e) \ 

" ^ /3ir 2 + a + r(M + e) D - \ 2 r 2 - X ± r + 9 + « - 7 r - ^r 2 ) 

(25) 

where D is the dilatation operator D — r J;. Transforming the system according to 

trW 0=o , t/=( C0SCJ " sin ^ (26) 

V 9 / l sinw cosu; / v ' 



with i? = y7 2 + u; = (1/2) arctan(/9i/7i). we obtain the following conditions for 
polynomial solutions to exist: 

Xl = R , Xl = P_É£±j2H (27) 

ix 

with 

deg p = n , deg q = n — 1 . (28) 
The final equations for p, q then read 

(D + A 2 )p + (A x r + A )q =0 
(/J + C 2 r 2 + dr + C )ç + (D x r + A>)p =0 (29) 



where 



. 2 (£-M) 03 l7o + ei? - Mi?) p 1 e + ^ R 

Al =7l ^^ + 71 ftï? + 

2 ct {aR — j3iK) i a(3± — kR 



A o =7it^ + 7i ^ ■+ 

C 2 =-2 

Ci = -2 



' Ai? i? 

C 2 = -2i? 

tfiM + 7071 



C 



i? 

7i« + 9R 
R 



, 2 M-e eR-MR-Pn , 7o i?- / 3 1 £, 

Di =-^iïjr +l1 — m — + ^?~ } 

^0 = ^ + ^^ + í*^ (30) 

Eq. (29) with (28) leads to a set of 2n + 3 algebraic equations in 2n + 2 parameters 
(9,E and 2n coefficients for p and q); so it cannot be solved for genèric vàlues of 
the physical coupling constants a, fli, 70, 71, M. However, considering one of these 
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physical parameters as free (e.g. the parameter 70) we are lead to consistent equations 
which, in principle, possess solutions and fix one energy level e and the constant 70 
as functions of the parameters a, 0i, 71, M and of k. 

The analysis of terms of highest and lowest degrees in r in these linear equations 
lead respectively to the following condition for the eigenvectors 

p = x n + o(x n ~ 1 ) , q = - T ^-^ x n - 1 + {x n ' 2 ). (31) 

-Ai 



(32) 

The first of these equations can be used to determine e while the second allows one 
to determine the parameter 6. 

We have analysed the equations numerically for the case n — 1, k — 1 correspond- 
ing to the ground state of the equation. Fixing symbolically 71 = 3,/?i = 4 and M — 1 
we found that two possible vàlues occur for the parameter 70. These vàlues are real 
for a G [0.097, 0.92] and become complex outside that interval. The two branches of 
vàlues are presented on Fig. 2. 

The results obtained above suggest that the reduced spherically symmetric Hamil- 
tonian is equivalent to an operator preserving the vector space P(n — 1) © P{n). 
However, in spite of our efforts, we could not find a change of basis making this job 
and we believe that it is not. Nevertheless, if we combine the two equations in order 
to obtain decoupled second order equations for the two components of the spinor, we 
got linear operators of the form [11] 

d 2 d 
T = (x 2 + x x) — + (-x 2 (x + x )+2p(x + x ))— + ex(x + x ) + (b-c)x + bx (33) 
dx z dx 

where 

£ = e 2_ M 2_ /£ _ 7 _ 1 ^ b = 2£a + ( K - 1 ) ; c= + (k-7) (34) 

a e + M a 

A necessary condition for the equation TQ n (x) = to possess polynomial solution of 
degree n in x requires è — n, quantizing the possible vàlues of the energy. Once this 
is fixed, the operator can be set in the form 

T = ~ x2 ( x -^ -n) - x x(x-^ -n) + (b- c)x + S QE s = xT QES + S QE s (35) 

where Tq ES , Sq E s can be expressed in terms of the three bàsic generators = 
x{xj- — n), = (x4- — | ) , J~ = 4-. This provides an alternative demonstation 
that T can be expressed as an element of the envelopping àlgebra of osp(2,2), as 
pointed out recently in [7]. Once set in this form, the equation TQ n (x) = leads to 
a system of n + 2 linear equations in n parameters (n parameters in Q n , remember 
that the parameter 7 is fixed by 7 = \/ k 2 — a 2 ). 



8 



So naively, one would expect the equations to fix two conditions among the cou- 
pling constants of the model, contrasting with the counting of equations and param- 
eters done with the systems of first order equations. 

However, due to the fact that the second order equation is obtained from the first 
order ones, it turns out that the two extra conditions are indeed consistent with each 
other and there is effectively one extra relation. To our knowledge, this property is 
not at all apparent by just looking at the second order equation. This construction 
can nevertheless be used to extend the class of one dimensional quasi exactly solv- 
able operators away from the class of operators which are directly expressible in the 
envelopping àlgebra of sl(2,R). 

5 Outlook 

We have shown that algebraic solutions of the Dirac + Coulomb + confining potential 
equations (with both Dirac oscillator and normal harmònic potentials) can admit 
some explicit bound states. A generalisation of the Dirac oscillator has been obtained 
which can be completely solved algebraically (i.e. it is exactly solvable) and posseses 
a hidden àlgebra related to osp(2,2). 

The mixed case, with both oscillators and Coulomb-terms, leads generally to an 
overdetermined systems of equations but, allowing one of the coupling constants to 
be "free", leads to a suflicient number parameters and the system of equations can 
be solved consistently by means of algebraic techniques. We believe that our results 
could be generalized by including screened potentials and using the ideas of [5]. 
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Figure 1: The vàlues of a and of the energy are plotted as functions of M in the case 
k = 1/2. 
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Figure 2: The two possibles vàlues of the parameter 70 are plotted as functions of the 
Coulomb coupling constant a for the "ground state" solution corresponding to mixed 
oscillators. 
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